Second Order Quasi-Normal Mode of the Schwarzschild Black Hole 



Hiroyuki Nakano 1 and Kunihito Ioka 2 
1 Center for Computational Relativity and Gravitation, School of Mathematical Sciences, 
Rochester Institute of Technology, Rochester, New York 14623, USA 
2 Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
(Dated: February 1, 2008) 

We formulate and calculate the second order quasi-normal modes (QNMs) of a Schwarzschild black 
hole (BH). Gravitational wave (GW) from a distorted BH, so called ringdown, is well understood 
as QNMs in general relativity. Since QNMs from binary BH mergers will be detected with high 
signal-to-noise ratio by GW detectors, it is also possible to detect the second perturbative order 
of QNMs, generated by nonlinear gravitational interaction near the BH. In the BH perturbation 
approach, we derive the master Zerilli equation for the metric perturbation to second order and 
explicitly regularize it at the horizon and spatial infinity. We numerically solve the second order 
Zerilli equation by implementing the modified Leaver's continued fraction method. The second order 
QNM frequencies are found to be twice the first order ones, and the GW amplitude is up to ~ 10% 
■ that of the first order for the binary BH mergers. Since the second order QNMs always exist, we 

can use their detections (i) to test the nonlinearity of general relativity, in particular the no-hair 
theorem, (ii) to remove fake events in the data analysis of QNM GWs and (iii) to measure the 
distance to the BH. 
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I. INTRODUCTION 



Thanks to the recent technological advance, we have almost come to the stage that gravitational waves are de- 
tectable. In the 21 century, the observation of gravitational waves will be absolutely a new window to our universe 
£~ ' i and also provide a direct experimental test of general relativity. 

There are several on-going projects for the gravitational wave detection in the world [l], 0, H, E| which are 
ground-based detectors. Next-generation detectors, such as the Large-scale Cryogenic Gravitational wave Telescope 
(LCGT) @ in Japan, are also in progress. In addition, as for space-based interferometric detectors, LISA 6] is now 
on its R & D stage and DECIGO/BBO HQ is proposed as a future project. Since they are space-based observation, 
they will be free from the seismic noise and remarkably sensitive to the low frequencygravitational waves below 1 Hz. 

One of the most important gravitational wave sources is ringdown of black holes [9j. The black hole perturbation, 
such as in the late stage of a black hole formation, can be described by quasi-normal modes (QNMs) with complex 
frequencies. Thus the gravitational radiation is expected as a damped sinusoidal waveform. It is important to study 
QNMs because we can determine the mass and angular momentum of a spining black hole by observing the QNM 
frequencies, i.e., the normal- mode frequencies and damping rates. 

For ringdown searches using data of gravitational wave detectors, the matched filtering technique is useful since 
these waveform is well understood. In the paper jic| . a data analysis method to search for ringdowns have been 
discussed by using an efficient tiling method for ringdown filters [Tljj , and an application to the TAMA300 data has 
been reported. Accuracies in the waveform parameter estimations have been found that [accuracy of black hole mass] 
< 0.9% and [accuracy of Kerr parameter] < 24% for events with the signal-to- noise ratio (SNR) > 10 [121 ]. 

A promising source that excites QNMs is a merger of binary black holes. In these events, we may detect the QNMs 
with high SNR, e.g., the SNR ~ 10 5 for ~ 10 8 M© black hole mergers at ~ lGpc by LISA [HI, since a large fraction of 
energy (~ 1-5% x [total mass]) is emitted as gravitational waves of the QNMs. Recently, numerical simulations have 
succeeded in calculating the entire phase of BH mergers 0EEEH, an d found that the £ = 2, m = ±2 mode actually 
dominates, carrying away ~ 1-5% of the initial rest mass of the system 17]. The merger rate is also estimated to be 
large enough [18|, [l9( . 

Black holes deform appreciably in the merger so that the higher-order QNMs could be prominent. As an order of 
magnitude estimate, when the gravitational wave energy of ringdown is ~ 1% x M, i.e., 

U(1)12 

Egw ~ M~ ~ 1%xM ' ( L1 ) 

where ip^ denotes the first order gauge invariant waveform function, the dimensionless amplitude of the metric 
perturbation is 

V~ lo% - (L2) 
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[See the total radiation energy of the first perturbative order in Eq. (|7.3ip with u ~ M Hence, the second order 
amplitude goes up to \ip^ 1 ' /Al] 2 ~ 1%. This means that the second order amplitude becomes ~ 10% of the first order 
amplitude, which is detectable for high SNR events. Reduced to its simplest terms, the SNR of the second order 
QNMs is - 10 if the SNR is ~ 100 for the first order QNMs. 

The essence of higher order QNMs has already been discussed in "Mechanics" by Landau & Lifshitz [20( as an 
anharmonic oscillation. In general, an oscillation with small amplitude x is described by an equation, 



with a frequency w. This gives a solution, 



x + u/x = (1.3) 



x = a cos(ut + 4>) , (1.4) 



where a and <f> denote integration constants. Including the second perturbative order with respect to the amplitude, 
the equation has a correction, 

x + uj 2 x = — ax 2 , (1-5) 

where a is a constant, and so the solution is x — acos(u;t + cf>) + x^ 2 \ where 

/,■, aa 2 aa 2 , ,„ . 

x {2) = -— - + — — cos2uitoc a 2 . (1.6) 
2w z bio 1 

Here the first term and second term in the right hand side of the above equation have a frequency (lo — lo) and (u> + u>), 
respectively. The important point is that the second order oscillation always exists and is completely determined by 
the first order one. We also note that the frequency does not have a correction at this order. 

In a recent paper [2lJ, we have proposed that higher perturbative order of QNMs, generated by nonlinear gravi- 
tational interaction near the BH, are detectable and worth searching in observations and simulations of binary BH 
mergers. In this paper, we give full details to calculate the second order QNMs of a black hole. Since the previous paper 
just put forward an order-of-magnitude estimate, we also perform numerical calculations of the second order QNMs. 
Considering a Schwarzschild black hole as a first step, we show that the second order QNMs appear at the frequ ency 
Ci/ 2 ) = 2a/ 1 ) and numerically determine their amplitude by the modified Leaver's continued fraction method [221 ]. 
Since the second order QNMs always exist, their detections can be used as a new test of general relativity. Although 
this paper and the previous paper [2l[ have been the first to study second-order QNMs, the second-order analysis of 
general relativity is pioneered by Tomit a [23l. a nd the I = 2, m = case is studied by Gleiser et al. [2~il. l25l 12a, l2l| . 
It is also extended to cosmology [H, [H, l30l. l3l| . 

This paper is organized as follows. In Sec. [ill we consider the second order metric perturbation and equations to 
be satisfied, i.e., the perturbed Einstein equation. Our strategy to solve this equation will be given in this section. In 
Sec. lIIIl we summarize the tensor harmonics expansion of the first and second order perturbation. In Sec. IIV1 we review 
the Regge- Wheeler- Zerilli formalism [32|, [33[ in the black hole perturbation approach. Here, the first order QNMs are 
also discussed for later use. In Sec.fVl we derive the second order Zerilli equation with a source term. The source term 
consists of quadratic terms of the first order wave-function. In Sec. IVI1 we regularize the second order source term 
so that it is regular at the horizon and spatial infinity. In Sec. I VIII we discuss how to extract physical information 
from the second order Zerilli function when we obtain this function. In practice, the gauge transformation from the 
Regge- Wheeler gauge to an asymptotic flat gauge is considered. In Sec. IVIII1 we numerically solve the second order 
Zerilli equation and calculate the QNM amplitude by implementing a modified Leaver's continued fraction method. 
In Sec. IIX[ we summarize this paper and discuss some remaining problems. Some discussions on the first QNMs 
are given in Appendix [A] We clarify the complex nature of QNMs and the relation between m and —m modes in 
the spherical harmonics expansion in this appendix. In this paper, we use units in which c = G = 1 and follow the 
conventions of Misner et al. [3J] with the signature — H — h + for the metric. 

II. SECOND ORDER METRIC PERTURBATION 

In the black hole perturbation approach, we consider second order metric perturbations, 

with an expansion parameter A which we can identify as a first order amplitude. Here, superscripts (i) (i = 1, 2) 
denote the perturbative order, i.e., h^} and hffi are called the first and second order metric perturbations, respectively, 
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and is the background metric. In this paper, we consider the Schwarzschild metric as the background and use the 
usual Schwarzschild coordinates, 



g^dx^dx" = -f(r)dt z + f(r)- L dr 
2 M 

/(r) = l-^. (2.1) 
r 

Thus, the Greek superscript and subscript indices denote {t, r, 6, <j)}. In the perturbation calculation, we raise and 
lower all tensor indices with this background metric. 

The Einstein tensor G^ v up to the second order is formally derived as 

GWBW = gW^VgW^I + g^^UWi+oi^)) 3 , j^w (^) 2 ), (2.2) 

where we have omitted the spacetime indices /i and v of the metric perturbations, h^v and hpv , and ignored the third 
perturbative order, i.e., 0(A 3 ) in terms of the expansion parameter. Gj}} is well known as the linearized Einstein 
equation, 

G^[H] = —-Hfiu-a'" + F^. u ) - R a ^ v H afi - -H.py, - -g^u{F\' x - H.\ X ) ; (2.3) 
F — H ;a 

Here, H^v denotes h$ or h$, and semicolon ";" in the index indicates the covariant derivative with respect to the 

(2) 

background metric. G\\v consists of quadratic terms in the first order perturbation, 

G$[hU,hW] = R^[h^,h^}-^R(V\h(% (2.4) 
rt®\hW - -h {1) hP-W +-h^ al3 (h {1) +h (1) -2h {1) ) 

2^ >p 2 a {f 1 ! u ) fjffia/ 2 m q ;p v 2 M a ;p v 

Since we consider the vacuum Einstein equation now, we may solve the following equation for the first perturbative 
order, 

= 0. (2.5) 

For the second perturbative order, once the first order metric perturbation h^ 1 ' is obtained, we may solve the equation 
with a source term which can be considered as an effective energy momentum tensor. 

G$[hW] = -G$[hV,h^]. (2.6) 

Thus, expanding the Einstein's vacuum equation, we can obtain basic equations order by order (35l. l36l. [3^| . 

In the following section, we consider the equations in Eqs. (|2.5p and (|2.6p by using the tensor harmonics expansion 
which is summarized in the next section. Then, the Regge- Wheeler- Zerilli formalism HH is used for coefficients 
of this expansion. 

In this formalism, the first order QNMs of a black hole is derived from the first order Einstein equation (|2.5p . We 
discuss the first order QNMs and summarize some formulae in the Regge- Wheeler- Zerilli formalism in Sec. IIV1 (See 
also Appendix [XI ) 

For the second order Einstein equation (|2.6[) , we also use the same Regge- Wheeler- Zerilli formalism in Sec. El while 
there are some differences. The second order equation has a source term which is written by quadratic terms of the 
first order wave-function. This source term which is shown in Eq. (|5.4p , does not behave well at the boundaries. 
The reason of the behavior simply comes from the gauge choice and this has no physical meaning. Hence, we can 
regularize this behavior in Sec. IVI1 Based on this regularized source in Eq. (|6.5p . we numerically compute the second 
order wave-function in Sec. I Villi 

In order to extract physical information from the first and second order wave- functions, we must consider the metric 
perturbations under an asymptotic flat gauge condition. The first and second order metric perturbations are obtained 
under the Regge- Wheeler gauge at first, but this gauge is not an asymptotic flat one. It is necessary to formulate the 
gauge transformation for both the first and second perturbative order. This is done in Sec. IVIII 
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III. TENSOR HARMONICS EXPANSION 



Since the background spacetime has the spherical symmetry, all perturbative quantities can be expanded by the 
spherical harmonics Y£ m (Q,<f)) and its angular derivative. In this paper, we consider the following tensor harmonics. 
Almost of all is the same as that of Zerilli's paper [33|. There are some differences in the notation, therefore we 
summarize them in this section. 

For the first and second order metric perturbations, we expand hfy (i = 1, 2) by tensor harmonics, 
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(3.1) 



where we should be careful not to confuse (t, r) in the above equation with the perturbed Einstein tensor G^l 
in the previous section. Here, ao£m, a £m , ■ ■ ■ are constructed by the spherical harmonics and its derivative. The ten 
tensor harmonics are defined as the following. 
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(3.9) 



(3.10) 



(3.11) 



y Sym - sin 2 9W em J 



Here the Sym denotes components derived from the symmetry of the tensors, and the angular functions Xi m and 
Wim are given by 



X 



Cm 



n d ( d 



1 d 2 



Ye m . 



(3.12) 
(3.13) 



From the above definition, the tensor harmonics can be further classified into even (or polar) and odd (or axial) 
parities. Even parity modes are defined by the parity (— 1) under the transformation (9,<fi) — ► (tt — 9, <fi + tt), 
while odd parity modes are by the parity (— l) e+1 . Thus, the seven coefficents of the metric perturbation, 

■^o'im' Hilmi ^2tmi ^vtm ^ltm' an< ^ ^tm are cane( i as tne even parity part, and the three coefficients, 

"O Cm 



i ^xtm ano - ^-2£m are ^ ne 00 -d parity part. 



On the other hand, we consider the right hand side of Eq. (|2.6p as the energy momentum tensor, 

Then we define the following tensor harmonics expansion, 

7"(2) _ ^ [A)^ m (t,r)a ^ m + .4u 7n (i,r)ai& n +^ m (i,r)a to + B £ m (t, r)b oem + Be m (t,r)b er , 



(3.14) 



Cm 



+Qoe m (t,r)c oem + Qem(t,r)ce m +T>i m (t,r)di m + Gi m (t,r)gt m + Tt m (t, r)ft m ] 



(3.15) 



Here the seven coefficents, Ao£ m , An m , Ai m , Boim, &im, Qim and Tn m are the even parity part, and the three 
coefficients, Qoe m , Qim and T>i m are the odd parity part. 

By using the orthogonality of the above tensor harmonics, we can derive the coefficient of the tensor harmonics 
expansion. For example, the coefficient of the energy momentum tensor is calculated by 



A £m(t, r) 



T^-a* oim dQ 



C M Q Cjy/3/T-(2) * 7Q 



(3.16) 



where * denotes the complex conjugate, d£l = sin 



and <5 MQ has the component, diag(l, 1, 1/r 2 , l/(r 2 sin 2 9)). 



IV. FIRST ORDER ZERILLI EQUATION 

We use the Regge-Wheeler-Zerilli formalism [32|, |33[ for the first order metric perturbation in the Schwarzschild 
spacetime. There are some reviews about this formalism in (38j . Separating angular variables with tensor harmonics 
of indices (£, m) as mentioned in the above section, the equations are divided into the even and odd parity parts. In 
this formalism, the master equation for the odd or even parity part arises as the Regge- Wheeler or Zerilli equation, 
respectively. 
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Here, we consider only the even parity mode in the first order calculation. The reason is the following. In the case 
of a head-on collision, we have already known that the odd parity perturbation does not arise due to symmetry. When 
we treat black hole binaries, the radial motion always exists in order to merge through the potential barrier of the 
system. Actually the numerical simulations show that the even parity mode dominates. Thus, the above assumption 
for the first perturbative order is adequate as a first step, though it may not be the best approximation. Furthermore, 
we should note that it is also sufficient to discuss the even parity part for the second order calculation under the above 
assumption. 

In order to discuss the metric perturbations, it is necessary to fix the gauge. For the first order metric perturbation, 
we impose the Regge- Wheeler (RW) gauge conditions, the vanishing of some coefficients of the first order metric 
perturbation: 

, (e)(l)RW _ , (e)(l)RW _ r {l)RW _ „ 
n 0£m —'him —^Irn ~ u ' 

Here, the suffix RW stands for the RW gauge. It is noted that in the RW gauge the gauge freedom is completely 
fixed. Although there are seven equations for the even parity part, introducing the following wave-function, 



A+l 



K. 



(1)RW 
tin 



(t,r) + 



r-2M 
A r + 3 A/ 



H. 



(1)RW, 



d 



■2 'in (*' r ) T Q r K tm 



(1)RW 



A = 



(£-!)(£ + 2) 



(4.1) 



we can reduce the seven equations to a single equation for the function V'im- This function ip^ obeys the following 
Zerilli equation. 



9 2 92 T7 t N 



V>i2(*,r) = 0; 



V z {r) 



2M\ 2 A 2 (A + l)r 3 + 6A 2 Mr 2 + 18AM 2 r + 18Af 3 
~) r 3 (Xr + 3M) 2 



(4.2) 



where 



r. = r + 2Mln(^_-l 



All the first order metric perturbations can be reconstructed from the wave-function V4rn' The detailed method of 
the first order 1 — 2 metric reconstruction under the RW gauge condition is summarized in the end of this section. 
If f/^m i s Fourier analyzed, 



(4.3) 



the Zerilli equation gives a one-dimensional scattering problem with a potential, 



drl 



Vz(r) 



(r) = 0. 



(4.4) 



Then, the QNMs are obtained by imposing the boundary conditions with purely ingoing waves, 



ib^ (r)e- luJt 

at the horizon of a black hole and purely outgoing waves 



-iuj{t-\-r* ) 



at infinity. Such boundary conditions are satisfied at discrete QNM frequencies. These frequencies are complex with 
the real part representing the actual frequency of the oscillation, i.e., the normal mode frequency, and the imaginary 
part representing the damping. There is an infinite number of QNMs for each harmonic index (I, to) which are labeled 
by n. Thus, the QNM frequency has three indices, (£, m, n). We note that the QNM frequencies have some symmetry 
shown in Eq. (|A"3l) . 
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In the following, we consider only the I = 2, m = ±2 modes as the first order perturbations. This is because these 
modes dominate for binary black hole mergers. Although we can derive all QNMs without fixing the m mode in the 
case of a Schwarzschild black hole, we need to specify m mode when we consider the second order perturbations. (See 
Appendix [A]) 

Furthermore, in this paper, we mainly concentrate on the most long-lived QNM in the first perturbative order. 
This mode is characterized by the fundamental (n — 0) QNM frequency. Just for comparison, we will calculate the 
second order QNM in the cases that only the n — 1 or 2 QNM is excited in the first perturbative order. Although the 
QNM frequencies for £ = 2 are given in [22j , we calculate them with higher-precision in Table U that is necessary to 
calculate the second perturbative order. 

For later use, we show that V4— 2 mus t be the complex conjugate of 4 22 m or der to assure that the metric 
perturbation has a real value. For example, the first order metric component h$ RW is written by 



Kt 



(1)RW _ 



fir) (H& 



f(r) (h^ 2 RW ^t)Y 22 {9^) + H$* 2 W (t,r) Y 2 _ 2 (8, 0)) 

' [V$ (*, r)} Y 22 (9, 0) + ff^T [4% (t, r)] F 2 _ 2 (6, <f>)) 
Ma \t,r)Y 22 {6,cj>) + 4%(t,r)Y 2 - 2 (0, 0)]) , (4.5) 



f(r) (H& 



where we have used Eqs. (|3.ip and (13. 2|) in the first line, and H^ 22 W — H^ 2 - 2 W is a real differential operator with 
respect to t and r defined by Eqs. (14. 8p . By using 



we can show 



Y 2 . 2 (9,cf>) = i\/-sin 2 #e- 2 ^, (4.6) 



~ 4e- , ^ , ('- r *) forr^oo, (4.7) 

where the second line represents the most dominant QNM (see Appendix [XJ and A is some complex number. Berti 
et al. 39j have discussed that the quasi- normal waveform for each (£, m) mode has four degrees of freedom, i.e., two 
amplitudes and two phases. In this paper, we consider only the even parity mode. This means that we need to 
consider two degrees of freedom in our calculation, i.e., the complex number A. 

From these wave-functions, the reconstruction of the first order metric perturbation under the RW gauge is given 
as the following, 

H^ 2 W (t,r) = H$™[^Ut,r)] 

2 r 2 - 2 r M + 3 M 2 d ,m , . 4r 3 + 4r 2 M + 6 rM 2 + 3 M 3 , m , , 
^2±2 (*> r ) ~ 3 777 77772 V>2±2 (t, r ) 



r(2r + 3M) dr^ ±2VW r 2 {2r + 3M) 

+ ( r -2M)-^4% (*,r) , 
Hi^(t,r) = H^ 2 W [^2 2 (t,r)] 

f 2 f) 

-4±2 (*,r) + r^-^ (t,r) 



2r 2 -6rM -3M 2 d , (i) °~ 
(2r + 3M) (r-2M) dt^ 2±2 {t ' T > + J 'dtd~r" p2±2 ' 



= H^ 2 W (t,r), 
-2M d 



r ~* M ° (tr) + 6 r 2 +rM + M 2 (1) 

In the above equations, we can derive the first order metric perturbation only from the wave-functions in simple 
differential forms. This is because there is no source term in the vacuum Einstein equation. On the contrary, when 
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we construct a second order metric perturbation, we need both the second order wave-functions and the second order 
source which have quadratic terms of the first order perturbation. 

We use the above metric perturbation under the RW gauge to derive the source term of the second order Zerilli 
equation. But, when we extract physical information, it is necessary to consider a gauge transformation from the RW 
gauge to an asymptotic fiat (AF) gauge. This will be discussed in Sec. IVIH 



V. SECOND ORDER ZERILLI EQUATION 



For the second order perturbations, we also separate angular variables by tensor harmonics and choose the RW 
gauge condition. First of all, we have considered the first order metric perturbation only for the even parity part. 
The second order source term, which is derived from G^l\h^ , h^] in Eq. (|2.6p . has also the even parity. Therefore, 
we may discuss the second order metric perturbation only for the even parity part, i.e., the Zerilli equation. 

Here if the dominant first order perturbations are the £ = 2, to = ±2 even parity part, we are sufficient to consider 
the £ = 4, to = ±4 even parity mode for the second order perturbations as shown below. The second order source term 
arises from G^j[/i (:L) ,/i (1) ]. This contains products, {I = 2, to = 2) x (£ = 2, to = 2), {£ = 2, to = -2) x (£ = 2, to = -2), 
(£ = 2, m = —2) x (Jt = 2, m = 2) and (£ = 2, to = 2) x (£ = 2, to = —2) of the first order metric perturbation. Then 
the source term has to = and ±4 modes as a result of the product of the spherical harmonics. We can find that 
the to = source term does not oscillate as a function of time, because of the symmetry of the QNM frequencies in 
Eq. (|A3|) . and hence this to = mode can not be observed as a QNM wave. Therefore, we consider only £ — 4, m = ±4 
modes in the following calculation. 

In [24], [2^, [27]], the second order source term of the £ = 2, m = mode which arises from the £ = 2, m = 
first order perturbations, has been discussed. These discussion is for example, for the treatment in the "close limit 
approximation" of the collision of black holes. Their situation is restricted to the axisymmetric case, i.e., the to = 
mode. Therefore, the second order calculation is also only for the to = mode. In our case, the harmonics of the 
second order metric perturbation changes in the to mode as well as in the £ mode as discussed above. 

Now, we introduce a function for the second perturbative order, 



2M 



3(3r + M) 



dK 



(2)BW 
4±4 



2M 



dt 



(5.1) 



where the functions K^2^ W arL d H'\1±a are coefficients in Eq. (|3.1j) . i.e., derived from the expansion of the second 
order metric perturbation by tensor harmonics as in the first order case. We note that the first-order counterpart 
exactly satisfies Xell = ^ipill/dt. Hence, the dimensions are i/}j^(t,r) ~ 0(M), Xim(t, r ) ~ 0(M°) and X£m(A r ) ~ 
0(M°). 

Using the function of Eq. (|5.1|) . the second order equations for the even parity mode is reduced to the Zerilli equation 
with a second order source term, 



r (2)RW 



dt 2 dr 2 



V z {r) 



*4±4(V' 



Si±i{t,r) . 



(5.2) 



where the potential Vz is the same function defined in Eq. 



with A = 9. The source term ^±4 is derived as 
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15(3 r + M) 9r"""' v " r ' 15 
4Vl0i7r (AM 2 + 7 Mr + 27 r 2 ) (r — 2 M) 



15 r (3 r + M) 2 



B oi±i (t,r) . 



(5.3) 



where the functions {$±±4 etc. are coefficients of the tensor harmonics expansion a£G^[h^\h^] in Eqs. (l3~14l) and 
p,15|) . Since these coefficients are written by quadratic terms of the first order metric perturbation, we can rewrite 
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the second order source term &t±4 in terms of V4±2 by using Eq. (|4.8[) as 

S i±i (t,r) = r — 1 - (24r 7 + 1120r 6 A/ + 1052r 5 A/ 2 - 798r 4 M 3 

' 42 v/il r 5 (3r + M) 2 (2r + 3M) 3 V 



-2586 r 3 Af 4 - 2396 r 2 AT 5 - 270 rM 6 - 9 Af 7 ) ( ^V4±2 (*,»")) J;V4±2 (*,»0 



+ = 7j ( 132 r 6 - 136 Mr 5 + 994 r 4 M 2 - 378 r 3 M 3 - 2306 r 2 M 4 

r 4 (3r + M) 2 (2r + 3M) 2 K 



-270 rM 5 - 9Af 



(78 r 4 - 128 r 3 M - 264 r 2 Af 2 - 100 rM 3 + 3 Af 4 ) / a 3 ( i) A d (1) 
r(2r + 3A/)(3r + A/) 2 (.-2A/) ^ 2±2 ( *' r) J d^ 2±2 ' 



(84 r 5 - 122 r 4 M - 281 r 3 Af 2 - 6 r 2 Af 3 + 237 rAf 4 + 33 M 5 ) / d 2 {1) \ d m 
r 2 (2r + 3M)(r _2M) 2 (3r + M) 2 ^ 2±2 ( '' r) J dt^ 2 



(66 r 4 — 106 r 3 Af — 220 r 2 Af 2 - 156 rM 3 — 45 A/ 4 ) / <9 3 (1) \ d (1) 



r (2 r + 3 A/) (3 r + M) 2 (r - 2 M) V ^ 2 9r ^ 2±2 ^ r) ) dt^ 2±2 ( *' " } 

-| j (336 r 8 + 3664 r 7 Af + 10144 r 6 M 2 + 15052 r 5 Af 3 



-13444 r 4 Af 4 + 7386 r 3 M 5 + 2648 r 2 A/ 6 + 270 rM 7 + 9 Af 8 ) 4± 2 r ) 7r^2±2 (*> r ) 



r 6 (3r + My (2r + 3Mf 

4- 9fU8 r- 2 A/f 6 -I- 970 rM 7 4- Q A/f 8 "\ i/i^ 1 ^ (i r\ 

at 

(r-2M)(3r + 2AJ) (1) \ _fl»_ (1) / fl» (1) \ (1) 

+2 r (3r + Af) 2 U 2 2±2 ( 7 ^ 2±2 ( j [dfidr^ 2 ^ ) drdt^ 2 ^^ 

1 5- (72 r 7 + 1240 r 6 A/ + 944 r 5 M 2 - 1284 r 4 Af 3 

r 5 (3r + Af) 2 (2r + 3Af) 3 v 

8 2 

-2910 r 3 Af 4 - 2396 r 2 Af 5 - 270 rM 6 - 9 A/ 7 ) v4± 2 (*, r ) g^W±2 (*> r ) 

(l32r* + 196, 4 Af + 174, 3 Af 2 + 114, 2 Af 3 + 28.M 4 -3^) < 1} (t>r) ^ ■ (5.4) 

r 2 (2r + 3Af) 2 (3r + Af) 2 (r-2Af) 2±2 V 7 St 3 r2±2 v ; J 

The second order metric perturbations can be reconstructed from x\±4 under the RW gauge as 

a K (2)RW (tr) 30r 2 + 9rAf + 2Af 2 (2) r-2M 9 (2) 

4±4 (t ' r) ~ r 2 (3r + Af) ^±4^^+ r g^A&r) 

zu/9 inrrfr- 9 M) 4^5 in r (r - 2 M ) „ . , 
Ui 4 ±4 (t, r) + , , 77 ; g 4±4 (* , r) , 







+ 3 " 


-(2)_rw 

2 4±4 


(i.r) 


= K< 
drdt 


.(2)W 
14±4 


(t,r) 


„3r + 
" r-2 


n(2)flW 

04±4 


(t,r) 


tt{2)RW 
— rt 24±4 ' 



3 r + M ^ v ' ' 15 3 r + M 

(t, r) + 3 ^-xi±4 (*) r ) ' 3 3r + M ^rxt±A (*j r ) - — H— «7Tf S 04 ±4 (*, r) , 



.J 2 ) ( f r ) + r2 9 K {2)RW (t r) 



8 V5 , 
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nr*T i±4 (t,r) . (5.5) 



Unlike the first order metric reconstruction, we need the information derived from G^f/i*- 1 ^, h^], such as Ai4±4 and 
So4±4 in Eq. (EUD. 
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VI. REGULARIZATION OF SOURCE TERM 

Although we show the second order i = 4, m = ±4 source term S4,±4 in terms of V4±2 m -^q. (|5.4[) , the raw expression 
54±4 does not behave well at infinity. This is not suitable for numerical calculations. We can find 5*4±4 ~ O(r ) at 
infinity because at large r the first order wave-function behaves as 



4%(t,r) = -F>'(t - u) + -Fj(t - r.) + [F,(* - r*) - MF'{t - r.)] + 0(r^) , 
3 r H 



(6.1) 



(2) 

where _Fj is some function of (t — r*) and FUx) denotes dFj(x)/dx. In order to obtain a finite solution X4±4i we 
have to make at least the second order source term ~ 0(r~ 2 ) by some regularization, which is the same order of 
the potential Vz ~ 0(r~ 2 ) in Eqs. (|5.2p and (|4.2p . On the other hand, the second order source behaves well at the 
horizon, i.e., ~ 0(r — 2Af), with 



$L&r) =F^(t + n) + 



lF H {t + r,) , 27F ff (t + r*) 



M 



56 



M 2 



(r - 2M) + 0[(r - 2M) 2 } , 



(6.2) 



where -Fy is some function of (t + r*) . 

We can regularize the second order source term by introducing the following regularized function, 



(2) 



C ^ it r) ^ (r ~ 2M)2 (±4,™ (t r)] JL^ (t 



•)• 



(6.3) 



The regularized function xi±I 69 satisfies the Zerilli equation (|5.2|) with a well-behaved source term, i.e. 



92 92 \t , ^ 



X { A±r 9 {t,r) = S i±i (t,r)- 



■w + m- Vz{r) 



(6.4) 



where S r ^ 4 ~ 0(r -2 ) at infinity and ~ 0(j~ — 2M) at the horizon. Thus we can remove an unphysical gauge- 
dependent divergence. 

It should be noted that such a regularization is not unique, and for example, we can replace d/dr with —d/dt 
in Eq. (|6.3[) . However the observed quantities do not depend on the choice of the regularization method. This is 

because the regularization is equivalent to adding quadratic terms in the first order gauge invariant function ip&±2 t° 

the second order gauge invariant function Ci±4; so that it preserves the gauge invariance [loT ]. 
The explicit expression for the regularized second order source term is given by 



S^^t, r) 



r — 2M V70 



1 



(72 r 8 + 3936 r 7 M + 2316 r 6 M 2 - 2030 r 5 M 3 



42 \ r 6 (3 r + M) 2 (2 r + 3 M) 3 

-7744 r 4 A/ 4 - 9512 r 3 Af 5 - 3540 r 2 M 6 - 1119 rM 7 - 144 M 8 ) (J^V4±2 <X r )) J^2±2 <X r ) 



1 



r 3 (3 r + M) 2 (2 r + 3 M) 2 (r - 2 M) 2 



(24 r 7 + 344 r 6 M - 872 r 5 M 2 - 771 r 4 M 3 



-120 r 3 A/ 4 + 77 r 2 M 5 - 237 rM 6 - ASM 7 ) ( ^-^±2 ) ^V4±2 (*,»") 



<% 2 

66 r 4 - 106 r 3 Af - 220 r 2 AT 2 - 156 rM 3 - 45 M 4 



',(1) 



dt 



(. 



d 3 



4% (t,r))f t i>& 2 (t,r) 



V drdt 2 



d 



(i) 



r (Af + 3 r) (2 r + 3 M) (r-2M) 
3 

'2160 r 9 + 11760 r 8 M + 30560 r 7 Af 2 + 41124 r 6 M 3 + 31596 r 5 M 4 



r 7 (2r + 3M) 4 (3r + M) 2 



+11630 r 4 Af 5 - 1296 r 3 M 6 - 4182 r 2 M 7 - 1341 rM 8 - 144 Af 9 ) t/4± 2 (*, r) ^^±2 (*» r ) 



r 5 (3r + Af) 2 (2r + 3M) 2 



228r' 



r°M - 370 r 5 M 2 + 142r 4 Af J - 384 r A M 



„3 ji#4 
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514 r 2 M 5 - 273 rM 6 - 48 M 7 ) (^42a (*,')) (*,r) 
(198 r 5 — 318 r 4 Af — 664 r 3 Af 2 — 458 r 2 Af 3 — 127 rAf 4 — 24 A/ 5 ) / <9 3 m \d m 



3r2(3r + M) 2 (2r + 3A/)(r-2M) ^ 2±2 ' 7 ^ 2±2 

' '216 r 8 + 4296 r 7 M + 1992 r 6 Af 2 - 3488 r 5 M 3 - 8716 r 4 Af 4 



r 6 (3r + Af) 2 (2r + 3Af) 3 

-9512 r 3 M 5 - 3540 r 2 Af 6 - 1119 rAf 7 - 144 Af 8 ) v4± 2 (*> ^g^2±2 (*. r ) 

5 g ( 252 r 6 + 636 r 5 Af + 674 r 4 Af 2 + 730 r 3 Af 3 + 524 r 2 Af 4 

r 3 (3 r + Af ) 2 (2 r + 3 Af ) 2 (r - 2 Af ) V 



+171 rM 6 + 24 M b ) ^ 2 (t,r) (*,r) + ^-^2 (^£2 (t,r) J ^^2 (t,r) j . (6.5) 

In the above equation, we note that the second order source term is quadratic in the first order wave-function ^±2 

and hence the second order QNMs have a frequency at lu^ — 2uj^ . Since the second order frequencies are different 
from the first order ones, we can in principle identify gravitational waves from the second order QNMs. 



VII. GAUGE TRANSFORMATION TO ASYMPTOTIC FLAT GAUGE 

In this section, in order to obtain physical information from the wave- functions, V4±2 an( i X4±I eS j we treat a gauge 
transformation in the first and second order calculation. To obtain the gravitational waveform, it is necessary to go 
to an asymptotic flat (AF) gauge from the RW gauge. This is because the Regge-Wheeler-Zerilli formalism that we 
have employed is under the RW gauge and this gauge is not asy mptotically flat. 

Here, we consider the following gauge transformation [H, l4l| . 

<w - Af = <w + GO + \ (* Q ) + (*")] , (7.1) 

where comma " ," in the index indicates the partial derivative with respect to the background coordinates, and 
anc [ £(2)m are generators of the first and second order gauge transformation, respectively. By using this gauge 
transformation, the metric perturbation changes as 



iW . 1,(1) _ 1,(1) _ r n (7 o\ 

a RWfiu ~ * n AFinv — a RW\iv ^P) V > ( ' • A ) 

1,(2) . .(2) _ 1,(2) _ 1 r „ , 1 r 2 „ r .(1) (7 o\ 



where C^i) means the Lie derivative in the £W direction. We use the following form of a generator for the gauge 
transformation, 

£«f = {V (i) (t,r)Y em (9,cl>), f(t,r)y fel (M), V 2 (i) (t, r)-^Y im (e, <j>), (t,r)-^^-Y lm (6^)} , (7.4) 

where we have considered only the even parity mode, and hence we have three degrees of gauge freedom for each 
order. (Note that the generator for the odd parity part is £Ma» = {0, 0, —Vgd^Ytm/ sin 6*, V3 svo.6 d$Yt m }.) In this 
paper, we consider only the t — 2, m — ±2 modes for the first order gauge transformation, and the £ = 4, m = ±4 
modes for the second order one. 

The gauge transformation of the metric perturbation is explicitly given as follows. For the first order metric 
perturbation, we derive 

H®g{t,r) = g&T (*■ r) + 2 §-V Q ^ (t, r) + 2 ^ ^ M) (t, r) , 
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H. 



K. 



(1)AF , 
2 2±2 

(1)AF 
2±2 



(t,r) 



^0 2±2 r ) 



h[ 2±2 (*! r ) 



G 



(1)^ 

2±2 



(t,r) = H^(t,r)-2^-V 1 ^(t, 



d 



(i) 



M 



r(r-2M) 



dr 

{ ^^V^it,r)-r^V^{t,r), 
r at 



Vi (1) (t,r) 



(t,r) = -2V 2 W(t,r) . 



(7.5) 



For the second order metric perturbation, we can calculate the gauge transformation straightforward, but obtain very 
long expressions. For example, they are written formally as 



(2)AF 



4±4 



K 



(2)AF 
4±4 



(t,r 
(t,r 



,l 14±4 V c ' ' 



where 5Hq 2 ^^, SK\± 4 , 5h ( il±l and 6Gf± 4 are defined by the tensor harmonics expansion of the last two terms in the 



,(2) 



,( 2 



^0 4±4 (*» r ) + a /0 ^ r ) + r(r _ 2M) 



-V^ 2 ' (t,r) + 6H$ ±4 (t,r) 



K 



(2)RW t 
4±4 



1. 



r 

,(2) 



{t,r)--y^ (t,r)+SK^ (t,r) 



= -^ 2 ( 2 )(t,r)+5Gf2 4 (t,r) , 



(e) 



,(2) 



(7.6) 
(7.7) 
(7.8) 
(7.9) 



right hand side of Eq. (73), i.e., (l/2)£ 2 (i)ff (o) _ £ 



. This includes only quadratic terms of the first order 



wave-function because hj^^ and are the first order quantities written by the first order wave-function ip. 
after we solve the gauge equations. 



(i) 

2±2 



About first order 



First, we consider the asymptotic behavior of the metric perturbation at large r in the RW gauge. Using Eq. 
and the asymptotic expansion in Eq. (|6.ip . the metric perturbation is given as follows, 



H, 



(1)RW 
2±2 



(t,r) 



H. 



22±2 
1 / rf 4 



H 



(1)RW 
1 2±2 



(t,r) 



3 v^ (T ^ r+ K^ ( ^) M+ ^ F/(Tr) 

(J^F I (T r )+ 1 4 (^ Fl (T r ))M+t (-^FriT^M* 

1 (_<£_ 
3 \dT r 
4 / d 4 



6 \dTr 3 'J 3 \dT r 4 

2 / d 4 \ Id 3 

- Fl {T r ))r--( -^Fj (T r ) j M — - ^F T (T r ) 



0(r~ 



3 \dT r 



K 



(1)RW 
2±2 



(t,r) 



3 \dT r 4 
1 d 3 

3 dT r 3 



Fi(T r ) )M 



dT r 



3 dT, 



•0(r" 



» + h (3 H " + sv F ' (T ->) + 0(r "» • (7 - 10) 

where we introduce T r = t — r*(r) for simplicity. Since we are interested in the asymptotic behavior now and the 
to — 2 and to = —2 modes have the same asymptotic behavior, we have used the same notation for the m = 2 and 
to = —2 modes. 

On the other hand, the metric perturbation in an AF gauge should behave as 

flS&Tfc r) = H®g(t, r) = h$m AF {t, r) = , 



J 02±2 

H^(t,r) = 0(r- 3 ) 



hmi AF {t,r) = 0(, 
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G&f (f,r) = Oir- 1 ). (7.11) 

This asymptotic behavior will be also used for the second order calculation. Then, we can find that the gauge 
transformation has the following form. 

^ M - -s (^< r '») r -5^ (r ' )+ K-535r ft < r '' + i GC^'W^" 2 '- 

tf, «- r) = 4^'< T ' ) 4?^' (T '» + M-^ M s|: F ' (7 '' ) -5 F ' (7 '' ) ) +0(O ' (7 ' 12) 

The above results are obtained iteratively for the large r expansion. Here, we note that the transverse-traceless tensor 
harmonics for the even parity part is fi m in Eqs. (|3 . 1 1) and (|3.11[) . Therefore the coefficient of the metric perturbation 
that is related to the gravitational wave is 

G hl AF - With E q s - COB and (I7T21) . we obtain 
G^2 F (t,r) = ll^ Fl (T r ) + 0(r- 2 ). (7.13) 
This can be shown in terms of V4±2 with Eq. (|6.ip as 

Ggg F {t,r) = ±4%(t,r) + 0(r- 2 ). (7.14) 



B. About second order 

In order to derive the gravitational wave amplitude for the second perturbative order, we also need to obtain the 
coefficient G)f±^ F under the AF gauge as in the first order case. From Eq. (|7.9|) . it is found that we must derive 8G^± A 
and V-^ 2 - 1 which is calculated by using the solution of V\" 2 \ In the following, we first obtain as Eq. (17. 20[) by 

solving the gauge equation in Eq. (17. 7[) . Then, is derived from Eq. (|T.8[) as in Eq. (|7.23p and 5G^± A is calculated 
as in Eq. ([7T25]) . 

(2) 

First we consider to obtain Vy in terms of the asymptotic wave-functions. The equation 

K<gf(t, r) = K { £* W (t, r) - V 2 ) (t, r) + SK^t, r) , (7.15) 

is used as shown in Eq. (|7.7[) . Here K*f±^ W is derived from xi±4i -4i4±4 and So4±4 in Eqs. (|5.5p . Since the wave- 
function xi±4 is given by xiil 69 and V 7 2±2 m Eq. (|6.3[) and the coefficients -4i4±4 and 23o4±4 are also given by 
t/4±2 with Eqs. (|3. 14[) . (|3.15[) and (|4.8|) . we may obtain Kf^ W in terms of wave-functions X4±T 9 an d ^2±2- The 

asymptotic expansion of i/4±2 ls given by Eq. (|6.ip while the asymptotic expansion of xli!* 69 is obtained from the 
Zerilli equation (|6.4|) as 

xi 2 2r(*,r) = FP (T r ) + 0(r-'), (7.16) 

(2) 

where Fj ' is some second order function. Here we consider only the leading behavior for large r, because this is 
sufficient in order to obtain the final gravitational waveform. If we consider the higher order calculation, it should 
note that there is a contribution from the second order source term of Eq. (|6.5p in the above equation. Then, we can 

calculate the asymptotic expansion of dK±±± W jdi in Eqs. (|5.5p as 
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M02y/^yy\dT r 5 y 7 dT r 3 v ; \dT, 

d 5 A d 2 „ . / d 4 „ A d 3 



dT r 

Integrating the above equation with respect to t, the second order metric component K±±^" W becomes 

K<g? W {t,r) = -jf>(T r ) 

70 / / d 4 ,„ A d 2 „ ,„ , / d 3 



34O20F I" VdT r 4 v r ')dT 2 K r> \dT, 3 



( d3 n ,m\\ ( d 



1-i 



\dT r 3Fl{Tr) ){dT? Fl{Tr)]}r]+()ir (7 "' LS) 

In Eq. (|7.15p . we also need the asymptotic expansion of SKW 4 , which is defined by the tensor harmonics expansion of 
^RWiiv m ^ C l- • Using the first order metric perturbation under the RW gauge in Eq. (|7.10p 
and the generator of the first order gauge transformation in Eq. (|7.4[) with Eq. (|7.12[) . 5K\±± is derived as 



From the above results and the AF gauge condition of in Eq. (|7.1ip . in order to remove the 0(r°) terms from 

Kfl* W and SKi% in Eq. (fTTTSjl . we obtain the pure second order gauge transformation V\ - ' as 

Vi (2) (t,r) = -r^ (2) (7V) 

70 A / d 4 „ \ d 2 „ ,„ , / d 3 ! 



r \ 9 -^TT-Pi — (T r ) +6 -_F z (T r ) 



13608^ ^ W?V 4 ) dT 2 \dT, 

+4 (^ F ' (T ' , )(^ F ' (T ' , ) ,u ) +0<rY <7 ' 20) 

Next, we derive V^ 2 -* from the AF gauge condition of h^l±l in Eq. (|7.11|1 . The second order gauge transformation 
of this component is given by 

h?m AF &r) = - 2( / 2M) ^ 2 > (t,r) - (t,r)+6h?m (t,r) , (7.21) 

as shown in Eq. (|7.8|) . Here we have used the fact that hfl±l RW = under the RW gauge. We may also derive 
Sh^l±l, which is defined by the tensor harmonics expansion of (l/2)C 2 (1) gffl — ^M^rwuu m ^1- Q7.3p . by Eqs. (|7. 10|) 
and with Eq. (fT~12"|) as 

^ = j^r ^ F , ( T rl+ (^ F , { T r) )^ + 0(,°). ,7.22, 

Then, V2' 2 ' is calculated from the above value and the result of in Eq. (|7.20p with the AF gauge condition in 
Eq. (j7TTj) as 

^ 2 )(t,r) = iif) (T r ) 

V70 1 /„„ / d 4 „ , m A d 2 . . on / d 3 



r(36 \—Fi{T r ) F r (T r ) +33 — ^ (T r ) 



136080? y VdT r 4 7 dT r 2 \dT, 

+ 4 (s£*''< r '>) {lf7 F - iT ' ) ) M ) +0 ^ ) 

^V 2 ,2) (t,r) + 0( r - 2 ). (7.23) 
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In the last line of the above equation, we have used the definition of T r = t — r*(r). 

At this stage, we can consider the metric perturbation related to the gravitational wave amplitude, i.e., G 4 ± 4 . 

(2) 

The gauge transformation of the components G\± A is given by 

G ( Hi F (t,r) = -V 2 ^(t,r)+6G { £ 4 , (7.24) 

as shown in Eq. (|7.9[) . Here <5Gjy| 4 , which is defined by the tensor harmonics expansion of (l/2)£| (1) g^ — C^d h^ W)ll , 
in Eq. (f773]) . is also obtained from Eqs. ([TTTOf and ([7^4]) with Eq. (|7TT2|) as 



Finally, using Eq. (17^51) for V 2 {2) in Eq. (17^1) . we obtain 

|o?ir«,r) = iff (IV) 

^;( 27 (^' (r ' , )^<-' +M &<- ) ) 2 
+4 (^ F '<^ ) )(^ F ' (I '' , ) M ) +0(r " ) ' (7 ' 26) 

By using V4±2 and X4±I eS j the gravitational waveform is obtained as 

+4M ^ t 4%(t,r^j ^vS 2 (*,r-)j + 0(r- 2 ) . (7.27) 
It should be noted that we can show that all metric components satisfy the asymptotic flat gauge condition. 

C. Power of Gravitational waves 

In this subsection, we summarize the power of gravitational waves which is derived from the gravitational waveform. 
The power P per solid angle is given by [12] as 

dP 1 / 1 / d , \ 2 1 / d , Id : '' 



dn ~ i^\ s ^e\dt hd *) + 4{m hee -^edi h **) /' (7 ' 28) 

where (• • •) means an averaging over a spacetime region which is sufficiently larger than the characteristic wavelength 
of gravitational waves. Therefore, for the £ — 2, m = ±2 modes of the first order metric perturbation, we have the 
following equation. 



+ ^G^ AF (t,r)W 22 (e,d>) + ^G^ F (t,r)W 2 - 2 (0,<f>) 



64tt 
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— C {1)AF (t 
g.' JT 2m K l i ' > 



+ ((|^(*.r) 



(x 22 (e, ( f > )) 2 + (w 22 (e,^)y 



l — n {x 2 _ 2 {e^)) 2 + (w 2 _ 2 (M)) a 



sin" 9 



' \x 2m (e,cf>)\ 2 + \w 2m (e,cf > )\ 2 



sin 

2 




1 (x 2m {9^)) 2 + (w 2m (e,4>)) 2 



sin"0 



(7.29) 



(1)AF 
22 



G 



(l)AF* 
2-2 



in terms of the coefficient of the tensor harmonics in Eq. (|3.ip . In the above equation, we use G 
and the angular functions Xg m and Wi m are given by Eqs. (|3.12[) . (|3.13p and (|4.6p . We note that it is meaningless 
to distinguish the power of m = ±2 modes since they appear as the cross term in the second equality of the above 
equation. It is also noted that the averaging ((dG 2 1 J n AF / dt) 2 } does not vanish in this calculation. Integrating the 
above equation with respect to the angular directions, we obtain 



p(i) 



8tt 



: E 

m=±2 



8tt 
_3_ 

4-7T 



E 

m=±2 



— G (1)AF (t r) 



(7.30) 



where we have used the result in Eq. (|7. 14[) in the second line and the fact that V4-2 = V*22 
line. The relation between m = 2 and m = —2 modes are summarized in the appendix. 



in Eq. (j4~7)) in the last 



When we set ipffi — A exp(— iw»= 2 m=2 {t — r *)), where we choose the amplitude A of ip 22 at the origin of t — r*, 
the total radiated energy is 



,(!) 



(1) 



= dtpw 



dt 



3 

4tt 
3 

47T 1 



c^ 22 



lAI 2 !^'! 2 / dt e 2 ^)(^r,) 



|2|,.,( 1 )|2_ 



1 



(7.31) 



For the second perturbative order, we also have the following formula for the power of gravitational waves from the 
result of the second order gravitational waveform in Eq. (|7.27[) , 



p(2) 



45 
8^ 

45 
4tt 



E 



m=±4 
xi4 )reS (^) + 




27 ^ (i, r) — ^ (t, r) + 24 ( (t, r) 



(7.32) 



Here, we have calculated the angular integration to derive the first line, and used the relation that the second order 
metric components of the m = 4 and —4 modes are complex conjugate, i.e., X44 T&9 = Xi-^ 9 * ■ This complex conjugate 
relation can be shown with Eqs. (|4.7[) and (|6.4p . We should note that when the first and second order perturbations 
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TABLE I: The first order QNM frequencies for the 1=2 mode. We have calculated them with higher-precision because it is 
necessary to calculate the second perturbative order. The QNM frequencies do not depend on m for a Schwarzschild black hole. 



n 




CV ,« 





0.37367168441804183579349200298 


-0.08896231568893569828046092718 


1 


0.34671099687916343971767535973 


-0.27391487529123481734956022214 


2 


0.30105345461236639380200360888 


-0.47827698322307180998418283072 



belong to the same harmonics {£, to) mode, there are the cross-term contributions even after the angular integration, 
such as G { H AF x Gfl AF * and G { ^ AF * x G^ AF from Eq. (OSl) . However, we do not have these terms in our case 
that G^ AF and G^ AF have different harmonics (£, to). 



VIII. NUMERICAL CALCULATION: LEAVER'S METHOD FOR SECOND-ORDER QNMS 



In this section, we consider a numerical method to obtain the particular solution X44 for the regularized 



second order Zerilli equation in Eq. 



Here we only consider the i — 4, m = 4 mode for the second order 



calculation because xi-I &9 mav be obtained by xi-I &Q = X44 reg * ■ I n order to solve the second-order Zerilli equation 
with several digits, we need to know many digits of the first order QNM frequencies. These are given in Table [U 

The following method is basically a modified version of the Leaver's continued fraction method [12]. We will finally 
calculate the second order amplitude at infinity and the horizon as the following forms, 



X44 reg (t, r = 00) = Ci L>ip£ (t, r = cx 
X $ re9 (t, r = 2M) = C H (t, r = 2M) 



(8.1) 
(8.2) 



where to denotes the first order QNM frequency ui^ an d we consider not only the fundamental n = mode, but also 
the overtone n = 1 or 2. And then, we may calculate the energy of the second order QNM as a function of the first 
order one, 



E&) (E^ 
= C e M\$uj 



where Eq. f7j32]) gives 



C E 



M 



20tt 



M 



Cr 



70 



15120F 



(51 - AiMuj) 



(8.3) 



(8.4) 



The coefficients, Cj, Ch and Ce, are summarized in Table HT1 
By Fourier transforming, 



we may solve the second order Zerilli equation with a source term, 

d 2 < -.x2 



where uj^ — 2u> (uj^ = 2u. 



0, 1 or 2) and the source term is given by 



(8.5) 



(8.6) 
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+3 
+ 



(276 r 7 + 476 r e M - 1470 r 5 Af 2 - 1389 r 4 M 3 - 816 r 3 M 4 - 800 r 2 M 5 - 555 rM 6 - 96 Af 7 ) oj 2 



9(r- 2M) 



r 7 (3r + Af) 2 (2r + 3M) 4 



r 3 (3 r - My (2 r + 3 Af ) (r — 2 M) 
2160 r 9 + 11760 r 8 Af + 30560 r 7 M 2 + 41124 r 6 Af 3 + 31596 r 5 Af 4 

(A(r)) 2 



T1630r 4 Af 5 - 1296 r 3 M 6 - 4182r 2 Af 7 - 1341 rAf 8 - 144 Af 9 



{r -2M) 2 M 2 u 2 



6(r- 2M) 



(144 r 8 + 4116 r 7 Af + 2154 r 6 M 2 



r 2 (3 r + Af ) 2 r 6 (3 r + Af) 2 (2 r + 3 Af ) 3 
2759 r 5 M 3 - 8230 r 4 Af 4 - 9512 r 3 Af 5 - 3540 r 2 M 6 - 1119 rAf 7 - 144 Af 8 



(r - 2Af) (7r + 4M)« 2 



3 (r — 2 M) 



r 5 {Zr + My (2r + 3M)' 



+142 r 4 M 3 - 384r 3 Af 4 - 514 r 2 M 5 - 273 rM 6 - 48 M 7 ) 



dr" 22 ^ 

(228 r 7 + 8 r 6 Af - 370 r 5 Af 2 



(8.7) 



-'22 



where = 

We solve the above equation under the boundary conditions with purely ingoing at the horizon and purely outgoing 
at infinity. This gives the two-point boundary value problem. However a simple numerical integration like a Runge- 
Kutta method does not work because the QNM frequency has a negative imaginary part and hence the QNMs diverge 
exponentially towards infinity and the horizon. We need unrealistic precision to resolve the unwanted solution, which is 
decreasing exponentially as we approach the boundaries. Furthermore, since the infinity is an irregular singularity f22j, 
the expansion of the basic equation around boundaries does not behave well, although this method was applied to 
the first-order [4311 . In the first-order case the most famous method to solve this problem is the Leaver's continued 
fraction method |22| . In the following we modify the Leaver's method and apply it to the second-order. 

In order to use the Leaver's method, we need to transform the Zerilli equation to the Regge- Wheeler one, and vice 
versa. This is possible by the Chandrasekhar transformations [H| for general (£, m) modes, 



V, (r) 



i 



A 2 (A + l) 2 + 4Af 2 w 



_ (2 6Af 2 (r-2Af) \ (1) 



2M^l(r). 



(8.9) 



We also have similar equations for X^Jj-jf ( r ) with replacing ui with a/ 2 ) — 2u> (0^44 = 2u>^; n = 0, 1 or 2). With the 
Chandrasekhar transformations, the Zerilli equation (|8 . 6[) becomes the Regge- Wheeler one, 



dr_ 

dr 2 



,(2) 



V RW {r) 



(2) reg,RW 

(2) 



(r) 
V RW {r) 



reg,RW i 



2Af\ (A + l)r-3M 
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r I 



(8.10) 



where the source term STw 2) (r) is also transformed by the Chandrasekhar transformations, but has a similar form 
asS£» ,(r)inEq. (EZD, " ' 



S 



reg.RW 

44a;( 2 ) 



(r) = (-0(^(r)) a + (-)^(r) f^f^ 



dr 



+ (••■) 



Here we can obtain the first order wave-function V^LM from ipzfa W (t) with Eq. ([8 

The Leaver's method provides the first order Regge- Wheeler function (r) m tne form, 



4 2 L W M = 2Af^(r)^a r , 



2M 



.11) 
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TABLE II: The coefficients for the amplitude and energy of the second order QNMs in Eqs. IpTI). ((8T2|) and (O) . The first 
order QNM has £ = 2, m = 2 with the overtone n. 



n CV 

0.221 - 0.489i 

1 0.397 - 0.384i 

2 0.459 - 0.226i 



Ch_ Ce 

-0.221 + 1.19i 15.0 
0.262 - 1.16i 12.7 
0.604 - 0.8Hz 8.97 



where 

^ V ; \2M J \2MJ 

p = -2iMw. (8.12) 

The coefficients a n are determined by three-term recurrence relations, 

a ai + (3aa = , 

a n a n+1 + (3 n a n + j n a n -i = 0, n = 1, 2, (8.13) 

where ao is arbitrary and a„, f3 n and 7 n are given by 

a n = n 2 + (2p + 2)n + 2p+l, 

(3 n = -(2n 2 + (8p + 2)?i + 8/9 2 + 4p + 2A-l), 

7„ = n 2 + 4pn + 4p 2 - 4 . (8.14) 
Similarly the second order solution can be written in the form, 

xZ 9 J W (r) = [A,(r)f £ ai? . (8.15) 

n=0 ^ ' 

(2) 

The coefficients a„ are determined by the following three-term recurrence relations, 

aoaf' + A)Oo 2) = b o > 
a n a% + n a,W + = b n , n = 1, 2, . . . , (8.16) 

where we replace p with 2p in a„, /3 n and 7„ in Eqs. (|8. 14|) and the source terms b n are determined by expanding 
S%j™ m Eq. (821 as 

- 3 l2 ^. /r-2Af\ n 



r-2M 

n=0 



We can first obtain the source terms b n algebraically, and then evaluate them numerically with the Maple calculator. 
Once we have numerical values of b n , we can solve the recurrence relations. The solution that satisfies the boundary 
condition at infinity is determined by adjusting a for which ^2 a n exists and is finite. 



IX. SUMMARY AND DISCUSSIONS 



In this paper, we have investigated the second order QNMs of a Schwarzschild black hole, as summarized in the 
following. 

1. We consider the £ = 2, m — ±2 even parity modes with the QNM frequencies £i4±2 f° r the first order perturba- 
tions because these modes dominate for binary black hole mergers. 
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2. We have considered the I = 4, m = ±4 even parity modes for the second perturbative order that are driven 
by the first order perturbations. We have employed the second order Regge- Wheeler- Zerilli formalism with the 
source term written by quadratic terms of the first order wave-function. 

3. We have regularized the source term because this does not behave well at the boundaries. 

4. Based on the regularized source, the second order wave-function have been obtained by using the modified 
Leaver's continued fraction method. 

5. We have explicitly derived the gauge transformation of the metric perturbation into an asymptotic flat gauge 
to extract physical information from the first and second order wave- functions. We have also formulated the 
power of gravitational waves for the first and second perturbative order. 

6. As the result, the second order QNM frequencies are found to be u4±4 = 2w2±2 an d the gravitational wave 
amplitude could go up to ~ 10%. This means that we can in principle detect and identify the higher order 
QNMs since their frequencies differ from any first order ones. The detectability of a two-mode ringdown wave 
has been discussed in [4a |. 

The detection of the second order QNMs would have the following advantages. 

1. The second order QNMs would provide a new test of general relativity, in particular of the no-hair theorem. 
According to the no-hair theorem, the astrophysical BHs are completely characterized by their mass and angular 
momentum. Then the mass and angular momentum derived from the second order QNMs should coincide with 
that from the first order ones. 

2. The first and second order QNM amplitudes include information about the total radiated GW energy E^ and 
E( 2 \ respectively. Using these ratio, this could provide distance indicators from a data analysis of only the 
QNM gravitational waves, since the observed GW amplitude is h ~ (E/M) 1 ' 2 (M/r). This is important for the 
case that the inspiralling phase is out of the detector frequency range. To be precise, since the observed GW 
amplitude also depends on the zenith angle 9, the parameters are degenerated and we can determine only the 
order-of-magnitude distances. A possible solution to this problem is presented below. 

3. We may also use the first-to-second amplitude ratio to reject fake events in the QNM search in which there are 
many fake events • In a method to remove fake events has been proposed by using the overtone QNM 
of the same (£, m) mode. The amplitude of the overtone QNM is determined by two factors, the 'excitation 
coefficient' of the QNM and the 'initial condition' of a perturbation. We can derive the excitation coefficient 
theoretically [46|, [43, [48|, [45^ but we must give the initial condition which depends on how QNMs are 
excited. Thus, it has been difficult to remove fake events by using the overtone QNM with an undetermined 
amplitude. On the other hand, the second order QNMs have the predictable amplitude except for the zenith 
angle dependence. This may ease the fake event rejection. 

Future problems include 

1. It is necessary to formulate the odd parity mode case in the Schwarzschild background. The odd parity mode 
appears when BHs have spin before mergers. 

2. We also need to discuss coupling contributions between different harmonics (£, m). 

3. As discussed in Sec. [V] the product of I = 2, m = ±2 even modes gives not only the m = ±4 modes, which is 
studied in this paper, but also the m = mode. Although the m = mode does not oscillate as a function of 
time, the amplitude would be comparable to that of the m = ±4 modes. Since the metric perturbation of the 
m = mode does not return to the initial value, this mode represents the gravitational memory effect [5ll l52l|. 
If the m = mode can be also detected, the amplitude ratio of the m = mode to the m = ±4 modes would 
determine the zenith angle to the observer since the zenith angle dependence is different. This could resolve the 
parameter degeneracy in the distance measurements with the second order QNMs discussed above. 

4. We have to extend the analysis to the Kerr BH case. When BHs have no spin before mergers, the final spin of 
the remnant BH is a ~ 0.7 QJ] and hence the Kerr effects may not be so large as inferred from the fact that 
the QNM frequencies shift by only a small factor. However, the final spin becomes large in the case of highly- 
spinning BH binaries |53fl. Since the master equation for Kerr BHs also has a source term that is quadratic in 
the first order function [54j| . we may expect similar results. 
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5. It is also important to discuss a mathematically rigorous definition of second order QNMs like the first order 
ones that use the Laplace transformation rather than the Fourier transformation l4q ] . 

6. The third order formulation is interesting. As suggested by [2(J, the QNM frequencies will also blueshift up to 
(-0W/M) 2 ~ 1% at this third order. 

7. In order to proof that the second-order QNMs actually exist and stand out of the GW tail, we have to find the 
second-order QNMs directly in the numerical simulations. Such simulations are challenging because the mesh 
size should be less than ~ 1% x M to resolve ~ 1% metric perturbations. There are detail analyses of nonlinear 
mode-coupling effects on a waveform by using the full nonlinear code [H, HE] > but the second-order QNMs have 
not been identified. 
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APPENDIX A: QNM FREQUENCIES 

In this appendix, we consider quasi-normal frequencies of a black hole in the first perturbative order in order to 
clarify their complex nature and relation between positive and negative m modes. First, we discuss the quasi-normal 
frequencies of a Kerr black hole, and then go back to the Schwarzschild case. 

With the spacetime symmetry, we may use the spheroidal (spherical in the Schwarzschild case) harmonics expansion 
which has the labels, I and m. There are infinite overtone frequencies of the QNM for each (£, m) mode. This is 
labeled by n. Even if we fix £, m and n, there exist two quasi-normal frequencies as we can see from Fig. 1 and 3 
in [221 ]. One has a positive real part Wt mn (+) while the other has a negative one Wg mn (_). 

For m > 0, the QNM frequency with a positive real part, cj£ m „( + ), represents more slowly damped modes than that 
with a negative real part, uji mn r_\, while for m < 0, uit mn r_\ represents more slowly damped modes than u>g mn ( + y 
Since we concentrate on the most dominant modes, we consider L0£ mn r + \ and o;^_ TO7l (_) (m > 0) in this paper. For 
example, although the integration in Eq. (|4.3j) picks up the two set of frequencies, uJe mn (+) and w^ rnra (_), we adopt 
the most slowly damped mode as 

iptm(t,r) = e^'-w'^^^W + e-^'-'Vta^,,., W (Al) 

f ^'VtaMwW. for m>0, 
I e - i ^o ( -)t^ ro ^ mo( _ ) ( r ) ! for m<0, 

where we consider the fundamental n — mode. We note that the QNM frequencies have the complex conjugate 
symmetry, 

Uimn(+) = — W|_ mn (_) , (A3) 

as discussed in Sec. 4 of [22] ]. 

Next, we discuss the Schwarzschild case. In this case, the m modes degenerate due to the spherical symmetry. 
Therefore, the label m have no meaning and we may write 

This equations means that the frequencies w^ m „( + ) and W£ m „(_) have the same damping time. The wave-function 
il>im(t,r) of the fundamental QNM is formally written by the same equation in Eq. (|A1[) . The two terms in the right 
hand side of Eq. (|A1|) have the same damping time in the Schwarzschild case, but we do not consider the term with 
W£ m „(_) for a positive m mode and LOi rnn (+) for a negative m mode, because these modes are not the more slowly 
damped mode in the Kerr case. Therefore, we have adopted the following wave-function for the first perturbative 
order, 

i>lm(t,r) = e _i ^-V £mw (i) (r) , for m>0, 

ili tm (t,r) = e^Jmty, ri>(r), for m<0, (A4) 
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where o;^ = UJ emn(+) an d we consider n = 0, 1 or 2 for comparison in this paper. These wave- functions satisfy 
Eq. (|4.7[) ensuring real metric components. 
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